April 2000 



Thermal Renormalons in Scalar Field Theory 

M. Loewe and C. Valenzuela 

Facultad de Fi'sica, Pontificia Universidad Catolica de Chile 
Casilla 306, Santiago 22, Chile 

Abstract 

In the frame of the scalar theory gcp^, we explore the occurrence of 
thermal renormalons, i. e. temperature dependent singularities in the 
Borel plane. The discussion of a particular renormalon type diagram 
at finite temperature, using Thermofield Dynamics, allows us to es- 
tablish that these singularities actually get a temperature dependence. 
This dependence appears in the residues of the poles, remaining their 
positions unchanged with temperature. 



1 Introduction 



During the last years, there have been an impressive amount of theoret- 
ical work on renormalons in different scenarios. For a recent review see 
One of the main motivations behind this effort concerns the non perturba- 
tive structure of Quantum Chromodynamics (QCD). On the other side, finite 
temperature effects have also called the attention of many authors , due to 
their crucial role in understanding thermal aspects of the hadron dynamics, 
whith special emphasis on the deconfining phase transition and the produc- 
tion of the quark gluon plasma 0. It seems natural, therefore, to start a 
systematic study of the occurrence and also the possible phenomenological 
role of thermal renormalons. Here we discuss, as a first step in this direction, 
the scalar theory g(f)'^. In spite of being a non realistic approach for phe- 
nomenological purposes, this analysis will allow us to gain a first impression 
on this kind of effects. 

The fundamental characteristic of a renormalon type diagram is an in- 
sertion of at least a chain of bubbles in a loop diagram, and behaves like 
k\ for large values of k, being k the number of bubbles in the chain. These 
diagrams, in the usual zero temperature situation, induce the existence of 
certain types of poles in the Borel plane. In this paper we explore in detail 
the extension to the finite temperature scenario of one kind of renormalon 
type Feynman diagrams, see Fig. 1, that contribute to the two-point func- 
tion. The set of diagrams we have chosen is an example of a renormalon type 
contribution (ultraviolet renormalons, UVR), and gives us a hint about the 
properties of the thermal Borel plane, i.e. positions and residues of poles as 
a function of temperature. 

The discussion at finite temperature has been done using the machinery 
of Thermo Field Dynamics (TFD) [Q. The main result from our analysis 
is that residues can get an explicit dependence on temperature, remaining, 
nevertheless, the position of the poles in the Borel plane unchanged. 

The plan of this paper is as follows. In section 2, the analysis at zero tem- 
perature of a particular set of Feynman diagrams, see Fig. 1, that contribute 
to the two-point function in the (^0^ theory, allows us to show the presence 
of a renormalon type singularity in the Borel plane. Section 3 is devoted 
to a brief discussion of how to handle the chain of bubbles shown in Fig. 2 
at finite temperature. This will be done in the deep euclidean region, i.e. 
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—p^ ^ m?, of the momentum p that circulates through the chain of bubbles. 
Using these results, in section 4 we calculate our renormalon type diagram 
(Fig. 1) at finite temperature and proceed to present our conclusions. 

2 The zero temperature renormalon 



In the present paper the diagram shown in Fig. 1 will be taken as a 
source for renormalons. Our goal is to explore the influence of temperature 
on this type of diagrams. Let us call by Rk{p) the diagram of order k, where 
k denotes the number of vertices. In order to establish our notation, we will 
first review the zero temperature calculation for the renormalon asociated to 
this diagram 0. 





Figure 1: Renormalon type contribution to the two-point function. 

If we denote by B{1) the one loop correction to the four point function, 
the so called "fish" diagram we have 

where the {—ig)^~'^ factor is due to the double counting of vertices in B{T). 
The relevant contribution to the integral comes from the deep euclidean 
region. Therefore, it is enough to approximate B{1) in this region. 

B{1) = tl^y f 



327r2 

(the argument of the logarithm is in mass units) 



{2ttY — rn? + ie {q + Z)^ — + ie 



log(-/^). (2) 



3 



In this way we have 



where a Wick rotation has been performed. 

This expression is uhraviolet divergent. Since g(j)^ is a renormahzable 
theory, we can concentrate on the finite part. In order to do this the propa- 
gator is expanded in powers of l/l^. The first two ultraviolet divergent terms 
are subtracted, and we keep only the first convergent term in the expansion. 
This procedure, which at the end induces one pole in the Borel plane, i.e. 
one renormalon, will be followed also in the next sections, where the finite 
temperature corrections will be computed. Due to this expansion, the de- 
pendence on p, the external momentum of the two point function, actually 
dissapears. So we find 

i?.oc-z(^)'/rf/i(log(/^))'^-^ (4) 

In the last expression we can see that the main contribution to Rk, at 
large values of k, comes from large values of I, that corresponds to the deep 
euclidean region, since we have done a Wick rotation. Introducing the new 
variable / = e*, we can see that Rk becomes proportional to the gamma 
function. 

Rk^-^{^^)'m■ (5) 

Let us remain breifiy the idea of the Borel transform. If, in general, we 
have a divergent series in terms of a certain expansion parameter a, of the 
form: 

oo 

D[a] = J2D^a^, (6) 

n=l 

then, one possibility to give a meaning to this series is to make use of the Borel 
transform method , according to which a new perturbative expansion, in a 
new expansion parameter b, is considered by dividing each coefficient of the 
previous series by n! in the following way 
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oo 



B[b] = Y.D. 



(7) 



n=0 



Formally, we see that the D[a] corresponds to the Laplace transform of 
B[b]. Using this idea in our case, we can introduce B[b] as the Borel transform 
of Rk according to 



The k\ behavior of Rk, for large values of k, is responsible for the ap- 
pearence of a pole in the real axis of the Borel plane, i.e. in the integration 
range of the Laplace transform. The meaning of this or other poles in the 
Borel plane corresponds to an ambiguity in the resummation of the series. 
Physically, this imply the existence of essential errors associated to pertur- 
batives expansions of physical magnitudes. 

3 Chain of bubbles at finite temperature 

In this section we revise the finite temperature calculation for the chain 
of bubbles shown in Fig. 2, which contribute to the four-point function. The 
original discussion can be found in |^. Later we will approximate our result 
for deep euclidean values of the momentum that goes through the chain. Let 
us denote by Ik the sum of all diagrams of the form shown in Fig. 2, with 
k bubbles and fixed external vertices of type one, according to the rules of 
TFD, including a global imaginary factor i. The sum is over all possible 
combinations of internal type of vertices, which, as we know, can be of first 
or second type. 

As it was shown in and after correcting some missprints, Ik can be 
expressed function of Ji according to 



m 




—I 



(8) 



oc 



l-6/327r2' 




(10) 



(9) 
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Figure 2: The chain of bubbles present in Fig. 1. 



where 



a(7) = 



A±iB 
9 



(11) 



In the previous expression a corresponds to the plus and 7 to the minus sign, 
respectively, and p is the momentum that circulates through the chain. A 
and B are given by 



A = Reli, 

pP\P0\ _ 1 

B ^ 



e^\po\ + 1 



I mil, 



(12) 
(13) 



So, for /i, the "fish" diagram at finite temperature including the global factor 
i, we have 



h = ii-igf 



2 J (271)^ [P - m'^ + ie 



+ 



27r5(P - 

pf^\lo\ - 1 



{I + py — m? + ie 



+ 



2t:5{(1 + pf - m'') 
e^Ko+po| _ 1 



(14) 



It is convenient at this point to give the explicit expressions for the real 
and imaginary parts of Ii^ and Ji^, zero and finite temperature parts of /i, 
respectively, such that /i = /i,, + Iig . 



Re I 



lo 



|4m2 — p"^ 



327r2 \ |p2 



log 



|4m^ _ p2| _|_ ^|p2| 



(15) 
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Re hp 
Im Ji„ 



9{p - 4m 



-9' 



Am? 



327r ' 
dll 



log 



{2p^Ef-{2l\p\+p 



2\2 



{2poEf-{2l\p\-p^y 



IQti\p\ Jo E{e(^^ - 1) 

"^'^ ^6{po + E, + E2) + 6{po-E,-E2 



8 J {2-kYEiE2 

+5{pQ + Ei- E2) + 5{pq -Ei + E2) 
/II 2 



X 



+ 



- l)(e/^^2 - 1)^ 



(16) 
(17) 



.(18) 



In the previous formulae Ei = y P + m^ and E2 = y {I + p)'^ + rn?, and 
the imaginary part has been obtained following |Q. Since a and 7 are con- 
jugated to each other, we can rewrite 



Reh = ga^cosikO), (19) 
Jm4 = ^aSin(/cg)(^ ^^|^^|_ J , (20) 
where we have introduced 



a = . (21) 
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9 = arctan(^^y (22) 

From the zero temperature calculation, we know that the renormalon type 
contribution of the diagram shown in Fig. 1 comes from the deep euclidean 
region, specifically for \f—p^ ~ e'^/^, where k is the number of bubbles in 
the chain, and p the momentum that circulates through the chain. In what 
follows, we will examine the chain of bubbles in the deep euclidean region at 
finite temperature, assuming a/— e^l'^ . Later, by replacing the exphcit 
expression for our chain of bubbles shown in Fig. 1, we will check if the 
previous condition, \f—p^ ~ e'^/^, is actually realized. 

First we calculate Imh^ and Re hp for T 7^ 0. From the four 5's that 
appear in Imhp, eq. |1^, it is easy to see that only those whose arguments 
include energy differences survive in the deep euclidean limit. So we have 
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By integrating and considering that e^'^^^ ^ 1, we have finally 

_ 2 -/3|p1/2 

Turning to the real part of Ji^, eq. note that the argument of the 
logarithm in the deep euclidean limit can be approximated in such a way 
that the integral can be written as 

where —p^ is sufficiently large so that we can reach the main contribution to 
the integral in the region where the logarithm can be expanded in powers of 
4/|p1/(-p2). In the previous expression, q denotes a certain bound for the 
integration in / such that for values of / bigger than g, the contribuition to 
the integral turns out to be negligible due to the exponential supression. At 
the end we can take q — * oo. The real part, then, can be written as a series 
in powers of 1/ \f--^ 

2 oo / ^n/2 

. (26) 

where the coefficients, which in our limit turn out to be essentially indepen- 
dent of the external momentum, are given by 

Using the previous results in eq. we can see that 

\kQ\ < 1, (28) 
where, once again, we have taken the assumption \/—p^ ~ e'''/^. Using this 



fact, we can approximate in eqs. [T^ and ^cos(/c^) ~ 1 and sm{k6) ^ 0. This 



means that Ik ~ ga ■ Since B in eq. 21 can be neglected then Ik ~ g{A/g) . 
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Therefore, using our expressions for the real part of Ji, eqs. and finally 

we get 



( 9 \^ 



(log( V))' + Mlog( + 



(29) 



This is the fundamental expression that we will use in the next section for 
the discussion of our diagram at finite temperature. Note that the leading 
or first term, which does not depend on temperature, is the same one we 
found in the zero temperature analysis. The second term is the first thermal 
contribution. The dots denote higher order terms in the expansion in powers 
of 

It is interesting to mention that we can find an approximated expression 
for the coefficients /„,(/?) in the low temperature region, where j3m ^ 1 

— 4"7r /-oo rjrf. „(n-l)/2 

^"^^^ 2 (n- 1)70 ^ ' 

which can be evaluated exactly, see [|I0| 



w) ^^i-X'^n"^) K.^m. (31) 



and where Ky are Bessel type functions 



4 Renormalons at finite temperature 



Using eq. ^ for the chain of bubbles, the renormalon diagram shown in 
Fig. 1 will be calculated in the frame of TFD. At finite temperature, Rk{T) 
denotes the sum of all diagrams of the shape shown in Fig. 1 with k vertices, 
being the external vertices of the first type. As it was the case in the chain of 
bubbles, the sum is over all possible combinations of internal type of vertices. 
In order to use our expression for J^, this sum must be performed before the 
integral over the internal momentum that circulates through the chain of 
bubbles. Here we will obtain an expression for Rk{T), finding the location 
of the induced poles in the Borel plane and the corresponding residues. 

The expression we have to deal with is 
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f dH ( i 27r5{ip + iy-m') 

J (27r)4\^(p + Z)2-m2 + ie e/'bo+'ol _ i 

The zero temperature contribution arises from the zero temperature part 
of the propagator times the zero temperature part of Ik-i- We would hke 
to discuss the following two cases: a) the product of the thermal part of the 
propagator times the zero temperature part of /fc-i, denoted by i?^, and b) 
the product of the usual zero temperature propagator times (including 
the zero temperature part of Ik-i), denoted by R^. 

Let us start with case a). We have 




J (27r)4 \ e^bo+'ol - 1 

The delta function excludes the deep euclidean region in the momentum 
I. In spite of this fact, we can exime the limit |[[ ^ m, after the integration 
in Z° has been done. It turns out that -Rfe(p) is proportional to 

The sum over k in the expression above is Borel summable and, therefore, 
it does not imply any renormalon. 

The second case, case b), corresponds to 



)iw^'^'°^^-''>n- '''' 



X 



1 \ 9^ 

-P 



+ 



(35) 



After the Wick rotation and substracting the divergents terms, exactly in 
the way as we did in section 2, we find that the last expression is proportional 
to 



327r2 




dt e 



-2t 



i2t)'-' + {k - l){2t) 



k-2-2t 



f2m 



, (36) 
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where the maxima of the first and second terms in the integral are reached at 
t = log(v— I^) ~ A:/2 and t ^ k/4, respectively. Note that the assumption we 
made in the previous section, \/—'p^ ~ e*^/^, has changed now in the second 
term, due to the factor 2 in the exponential. This fact does not affect our 
approximations. 

Doing the integral, we see that i?^ is proportional to 

) ^(fc)-2^/2(/5)(^) r(fc) + --.. (37) 

Wc sec that those terms proportional to the product of the tempera- 
ture dependent part of the propagator times the whole series of are Borel 
summable, since the leading term in the series of 7^, case a), is already Borel 
summable. On the contrary, the zero temperature part of the propagator 
times /fc gives us a series of terms that behave like k\ being, therefore, non 
Borel summable. We have calculated explicitly the first temperature depen- 
dent term of this series associated to a k\ behavior. 

The Borel transform, of the sum Yl^Rk, taking into account only 

those non-summable Borel terms, is given by: 



Rl oc 



Rk 



Lk-l 



oc — z 



6/327r^ 



+ 2/2(/3) 



6/647r^ 



+ 



(38) 



The first term in the last expression is the zero temperature rcnormalon 
we found in section 2. The second term corresponds to the first thermal 
singularity in the Borel plane, and the dots denote a whole series of poles. 
The structure of singularities in the Borel plane, associated to our renormalon 
type diagram, is shown in Fig. 3. The leading or first pole along the real 
axis is the zero temperature renormalon. The series of new thermal poles, 
i.e. thermal renormalons, to the right side of the first one, are subleading 
contributions to the Borel sum ambiguity. 

Summarizing, in this paper we have discussed thermal corrections to a 
particular renormalon type diagram in the theory g^^^. Prom this analysis, we 
show the existence of thermal singularities in the Borel plane associated to 
this particular diagram. The main properties of these thermal renormalons 
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A 



X X X — >f- 



T=0 



Figure 3: Singularities in the Borel plane. 

are the following: a) Their location in the Borel plane does not depend on 
temperature. They are situated at the points where the zero temperature 
renormalons are usually located, i.e. at n/Po, being n an integer and Po the 
first coefficient of the (3 function |TI|. b) Their residues, on the contrary, 



have an explicit dependence on temperature, through the factors /„(/?) that 
vanish when T — 0. 

The conclusions from the particular diagram we have discussed here, sug- 
gest us to conjecture that these will be actually general properties of the Borel 
plane associated to correlation functions at finite temperature: residues, in 
general, will depend on temperature being, however, the location of the poles 
temperature independent. Thermal renormalons in the scalar theory are al- 
ways subleading in the ambiguity of Borel sum, i.e. the leading pole is tem- 
perature independent. This fact is related to the ultraviolet character of the 
renormalons in this theory. Thermal corrections are associated to long dis- 
tance correlations in the system, whereas the leading renormalons come from 
the singular behavior at short distances. Something equivalent happens with 
the axial anomaly, where it has been shown that the anomaly itself does not 
depend on temperature |T^ . In QCD this situation will be probably different 



due to the existence of infrared renormalons. 
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